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Network  Identification  Performance: 

(1)  Simulations  for  the  Middle  East/North  Africa,  and 
(2)  Maximum  Likelihood  Estimates  of  Teleseismic  mb  for  the 
GSETT-3  Primary  Network 

Terrance  G.  Barker 


1.  Simulations  for  the  Middle  East/North  Africa 

1.1  Introduction 

We  have  simulated  the  detection  and  identification  performance  of  the  current  and 
proposed  IMS  seismic  networks  in  the  Middle  East/North  Africa.  Figure  1  shows  a  map 
of  Africa  and  western  Eurasia  with  the  current  IMS  alpha  stations  and  the  area  for  which 
the  simulations  were  done.  The  identification  performance  of  a  network  is  strongly 
dependent  upon  regional  source  and  propagation  variability.  However,  knowledge  of 
those  variations  allow  one  to  estimate  their  effects  and  to  know  that  in  some  areas  certain 
discriminants  can  work  and  in  other  areas  they  do  not.  Data  sources  included  published 
results  on  wave  propagation  in  this  region  and  incorporated  analyses  of  data  taken  by 
Vernon,  et  al,  (1996)  in  Saudi  Arabia. 

1.2  IGPP  Data  Set  and  Regional  Properties 

We  recently  received  data  collected  by  Vernon,  et  al.  (1996)  on  the  Arabian  Shield. 
The  data  were  provided  to  us  by  Frank  Vernon.  Using  the  larger  events  in  their  data  set 
for  which  there  were  magnitudes  reported  in  the  IMS  Reviewed  Event  Bulletin,  we 
extracted  parameters  which  we  used  to  make  simulations  for  this  region.  The  stations  in 
their  array  and  the  events  we  analyzed  are  plotted  in  Figure  2.  The  IMS  coverage 
currently  is  quite  sparse  in  this  region,  so  their  data  represents  an  important  increment  in 
knowledge  there.  In  addition,  the  data  set  allows  us  to  infer  noise  levels  directly  for 
stations  on  the  shield. 

Before  analyzing  specific  events  in  the  dataset,  we  examined  the  entire  origin  file 
to  infer  an  approximate  mb  threshold  for  the  region.  If  the  probability  of  detection  is  a 
Gaussian  distribution  P(mb;  p,o)  with  mean  p  and  sigma  o,  then  we  assume  the  observed 
distribution  follows 

Ninc  =  P(mb;fi,a)N0CC(mb), 

log(Nocc)  =  a-b\ogmb 

We  further  assume  that  b=l  and  then  fit  the  distribution  of  observed  magnitudes,  which 
were  from  IMS  bulletins.  Figure  3  shows  the  data  and  fit,  with  a=6.5,  p=4.5  and  0=0.3. 
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Figure  1.  The  current  IMS  network  and  the  area  for  which  the  simulations  were  done. 

Thus,  the  threshold  for  the  IMS  at  the  time  these  data  were  taken  (late  1995  to  early 
1996)  for  the  region  as  a  whole  was  over  mb=4.5.  By  running  simulations  for  that  network 
we  were  able  to  set  the  scale  factor  in  the  P  wave  attenuation  curves  (effectively  setting 
the  relationship  between  moment  and  mb). 
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Figure  2.  The  network  used  by  Vernon,  et  al  (1996)  (solid  triangles)  and  the  locations  of  earthquakes 
(open  circles)  analyzed  for  this  study. 

Vernon,  et  al,  show  that  events  in  this  region  are  very  different  in  terms  of  their 
magnitude  and  relative  amplitudes  of  regional  phases.  We  selected  events  to  represent  five 
distinct  areas,  as  suggested  by  Vernon,  et  al.,  and  measured  time  domain  and  spectral 
values.  In  order  to  source  excitation  levels,  we  calculated  log(Lg/P)  values  and  log 
amplitudes  corrected  to  mb=4  at  a  distance  of  1000  km,  relative  to  events  in  the  NE 
Arabian  Shield.  That  is, 

log  Acorr  =  log  Aobs(mb,  A)  +  mb  -  4  -  0.833  log(  A  / 1 000)  (1) 

It  would  be  preferable  to  relate  the  values  to  a  long  period  measure,  but  long  period 
spectral  values  (computed  by  us)  from  some  of  the  areas  were  poor  and/or  there  were  not 
Ms  values  from  each  of  the  regions.  It  would  also  be  preferable  to  calculate  sub-regional 
values  of  Q(f),  but  this  will  be  done  by  other  institutions.  Here,  we  assumed  that  values 
above  200  for  Lg  Q  on  the  shield,  near  200  off  the  shield  and  very  low  in  the  Red  Sea. 
Table  1  shows  the  observed  relative  amplitudes. 
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Figure  3:  Observed  distribution  of  events  with  mb(+)  and  fit  to  the  distribution(x). 


Table  1.  Observed  Relative  Amplitudes 


Source  region 

Lg/P 

Acorr  f°r  L9 

A  for  Pn 

NE  Shield 

0.3 

0. 

0 

S.  Red  Sea 

1 

-0.8 

-1.5 

Aqaba 

2 

0.5 

-1.1 

Zagros 

-0.3 

-1.5 

-0.9 

Arabian  Sea 

0.4 

-1 

-1 

The  study  region  was  made  into  a  grid  that  represents  differences  in  source  and 
propagation  properties  as  reported  by  Vernon,  et  al.  and  summarized  by  Sweeney  (1995). 
The  values  of  Table  2  and  the  other  studies  were  used  to  populate  the  grid,  which  is 
shown  in  Figure  4. 

For  each  of  the  grid  members,  we  used  relative  amplitudes  of  signals  from  the  sub- 
regions  to  set  the  regional  parameters.  In  Xnice,  the  spectral  amplitude  of  a  regional  phase 
at  distance  A  has  the  form 
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Figure  4.  Grid  of  source  and  propagation  parameters  used  in  the  simulations. 


Licind 
T  tectonic 

AF  ifrica 
R3  rad  saa 
AS  arablan  saa 
ASH  arablan  shlald 
Z  zaoros 

AQ  acaba 
LG  B  Ic  blk 


log  Aeq  =  (f)+S0-  0.8331og(  A  /  A0 )  -  0.4342y(A  -  A0 ), 

log  Aa  =  S„(f)  +  S0~  0.833  log( A  /  A0)  -  0.4342y(A  -  A0)  +  C0exp  +  C,exp  log(/),  (2) 

where  Seq(f)  and  Sa  (/)  are  the  earthquake  and  explosion  source  spectra,  and 

r  =  nf^Ug-'Q0-\ 

Here,  U ’  is  group  velocity,  Q0  is  the  quality  factor  at  1  Hz,  r|  is  the  slope  of  the  Q 
versus  frequency  relation,  S0  is  the  source  strength  for  unit  moment,  and  C  are  the 
explosion  relative  source  strength  coefficients  (see  Barker,  1996  for  details).  The 
earthquakes  generally  occur  in  a  narrow  distance  band  (800  to  1100  km)  from  the 
network,  so  it  was  not  possible  to  separate  source  and  propagation  factors  from  these 
data  alone.  We  relied  on  the  summary  in  Sweeney  (1995)  for  Q  values  and  used 
Baumgardt  (1996)  and  Vernon,  et  al.,  to  deduce  relative  source  strengths.  There  are  no 
explosions  in  the  data  set,  so  we  arbitrarily  set  C0exp  to  zero  for  Lg  and  to  one  for  Pg  and 
Pn.  Since  we  could  not  reliably  estimate  C,exp,  we  did  not  examine  discrimination 
performance  based  on  explosion  spectral  slope.  The  parameters  used  in  the  simulations 
discussed  here  are  summarized  in  Table  2. 
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Table  2.  Regional  Source  and  Propagation  Parameters 


africa 


Phase 

Qo 

il 

Co'*' 

c,cxp 

Ls 

3500. 

600. 

0.39 

-20.9 

0. 

-0.1 

Pg 

6000. 

600. 

0.39 

-21.2 

1.0 

-0.5 

Pit 

8500. 

600. 

0.6 

-21.2 

0.8 

-0.5 

aqaba 


Phase 

u. 

Qo 

T| 

S» 

C0‘xp 

C,exp 

Lg 

3500. 

200. 

0.39 

-20.4 

-0.1 

pg 

6000. 

200. 

0.39 

-22.4 

-21.9 

-0.5 

Pn 

8500. 

200. 

0.6 

-22.4 

-21.9 

-0.5 

Phase 

U, 

ara 

Qo 

bian_s« 

T1 

;a 

C0'xp 

C,exp 

Lg 

3500. 

200. 

0.39 

-21.9 

0.0 

-0.1 

Pg 

6000. 

200. 

0.39 

-22.3 

1.0 

-0.5 

Pn 

8500. 

200. 

0.6 

-22.3 

1.0 

-0.5 

arabian_shield 


Phase 

Us 

Qo 

11 

S0 

C0'xp 

C,cxp 

Ls 

3500. 

300. 

0.39 

-20.9 

0.2 

-0.1 

Pg 

6000. 

300. 

0.39 

-21.2 

1.4 

-0.5 

Pn 

8500. 

300. 

0.6 

-21.2 

0.8 

-0.5 

lgjblk 


Phase 

U, 

Qo 

11 

So 

C0cxp 

c,cxp 

Lg 

3500. 

5. 

0.45 

-21.7 

0.4 

-1.1 

Pg 

6000. 

10. 

0.68 

-22.7 

0.8 

-0.1 

Pn 

8500. 

10. 

0.8 

-22.7 

2.2 

-0.1 

red_sea 


Phase 

Uc 

Qo 

O 

So 

c0cxp 

C  cxp 
^1 

Ls 

3500. 

40. 

0.45 

-21.7 

0.0 

-1.1 

Pg 

6000. 

40. 

0.68 

-22.7 

1.0 

-0.1 

Pn 

8500. 

40. 

OO 

d 

-22.7 

1. 

-0.1 

tectonic 


Phase 

Us 

Qo 

11 

So 

c0exp 

C,exp 

Ls 

3500. 

671. 

0.43 

-21.2 

0.3 

-1.1 

Pg 

6000. 

1000. 

0.39 

-21.5 

1.4 

-0.1 

Pn 

8500. 

500. 

0.6 

-21.5 

0.8 

-0.1 

zagros 


Phase 

U, 

Qo 

11 

s. 

c0exp 

c,cxp 

Ls 

3500. 

200. 

0.39 

-22.6 

0.0 

-0.1 

Pg 

6000. 

200. 

0.39 

-22.3 

1.0 

-0.5 

Pn 

850071 

200. 

0.6 

-22.3 

1.0 

-0.5 
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1.3  Noise  Levels  of  Saudi  Stations 


For  the  simulations  used  here,  we  used  time  series  from  the  IGPP  data  set  to  find 
noise  levels  for  the  proposed  SAUD  station.  The  Saudi  stations  are  very  quiet,  as  can  be 
seen  from  Figure  5,  which  shows  representative  noise  spectra  for  five  stations  on  the 
same  day,  plus  the  low  noise  model  proposed  by  Peterson  (1993).  The  station  noise 
levels  are  comparable  to  the  low  noise  model,  with  the  noise  at  station  HALM 
significantly  lower. 


Figure  5.  Noise  spectra  for  five  Saudi  stations  on  the  same  day  (1995339),  plus  the  low  noise  model 
proposed  by  Peterson  (1993). 


1.4  Simulations  of  Detection  and  Identification  Performance 

The  IMS  network  proposed  by  the  Conference  on  Disarmament  Experts  working 
Group  in  August,  1995  is  more  dense  in  the  Middle  East  (Figure  6)  than  that  which  is 
currently  in  place.  In  the  following  we  compare  estimates  of  the  performance  of  the 
existing  network  with  the  proposed  net.  We  begin  with  the  detection  thresholds. 
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Figure  6.  Proposed  IMS  seismic  stations. 

Contours  of  mb  detection  thresholds  for  the  current  network  are  shown  in  Figure 
7.  By  detection,  we  mean  the  criteria  used  in  GSETT-3  and  currently  implemented  on  the 
IMS.  The  rules  are  based  on  a  weighted  sum  of  measurements  (J.  Carter,  Center  of 
Monitoring  Research,  personal  communication): 


m  ft  Jtq.e  dfiS  .nor»g.mb  Pthreeh 


30’  45'  60‘  75‘ 


30'  45*  BO'  75' 


Figure  7.  Contours  of  mb  detection  thresholds  for  the  current  network. 


Table  3.  Weights  for  Measurements  in  Determining  Detection  Status 


Primary  Time 

1.0 

Secondary  time 

0.7 

Array  azimuth 

0.4 

Array  slowness 

0.4 

Single  station  azimuth 

0.2 

Single  station  slowness 

0.2 

Her,  “primary”  refers  to  teleseismic  P  or  regional  Pg  or  Pn  and  secondary  refers  to 
teleseismic  pP  or  S  or  regional  Lg.  All  azimuths  are  determined  from  P,  Pg  or  Pn.  A  sum  is 
formed  from  the  number  of  measurements  which  exceed  the  signal-to-noise  ratio  times 
their  respective  weights  for  the  stations  in  the  network.  An  event  is  considered  detected  if 
the  sum  exceeds  3.55.  Thresholds  in  this  figure  and  in  the  remainder  of  this  report  are 
90%  confidence  values. 

The  thresholds  decrease  dramatically  towards  the  NE  as  coverage  improves.  For 
events  on  the  Arabian  peninsula,  detections  are  teleseismic  (no  stations  are  within  the 
regional  cutoff  distance  of  20  degrees).  The  detection  thresholds  for  the  proposed 
network  are  shown  in  Figure  8. 
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Figure  8.  Contours  of  nib  detection  thresholds  for  the  proposed  IMS  network. 


The  additional  stations  lower  the  detection  threshold  substantially.  Variability 
within  the  region  is  due  to  station  geometry  and  variations  in  source  strength  and 
propagation  parameters  (Tables  1  and  2,  Figure  4). 

We  turn  now  to  identification  performance.  We  note  first  that  teleseismic 
discriminants  are  ineffective  because  events  in  this  region  (and  the  simulation)  are 
typically  shallow  (<10  km).  We  did  not  model  the  Ms:mb  discriminant  because  we  need 
accurate  estimates  of  long-period  noise  for  the  networks,  which  we  currently  do  not  have. 
Thus,  for  events  occurring  on-shore,  we  must  rely  on  regional  discriminants  for 
identification.  Since  there  are  no  events  detected  (in  the  sense  of  the  rules  in  Table  3)  by 
the  current  IMS  network  in  this  region,  discrimination  cannot  be  done.  On  the  other  hand, 
the  proposed  network  is  capable  of  regional  discrimination.  We  consider  first  the  Lg/P 
discriminant.  An  event  is  considered  identified  as  an  earthquake  if  (1)  the  event  satisfies 
the  event  detection  criterion  (Table  1),  (2)  both  Lg  and  Pn  or  Pg  exceed  the  signal-to-noise 
ratio  at  least  1  station,  and  (3)  the  Lg/P  ratio  exceeds  a  specified  value.  Figure  9  shows  the 
contours  of  mb(Lg)  threshold. 
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Figure  9.  Contours  of  mb(Lg)  detection  thresholds  for  the  proposed  IMS  network. 

Note  that  to  the  north  of  SAUD,  the  mb(Lg)  threshold  increases  locally.  This  is 
due  primarily  to  the  reduced  excitation  of  Lg  (see  Table  1)  towards  the  Zagros 
propagation  grid  zones.  We  can  see  the  influence  of  the  excitation  levels  So  on  nib(Lg) 
from  Figure  10,  which  shows  contours  those  levels.  In  general,  where  the  excitation  is 
high,  the  threshold  is  low,  and  vice-versa.  In  addition  to  the  detection  of  Lg,  the  Lg/P  ratio 
must  exceed  a  specified  level.  We  show  in  Figure  10,  contours  of  observed  Lg/P,  which  are 
low  towards  the  Zagros  region,  and  high  in  the  Gulf  of  Aqaba.  We  would  thus  expect 
from  Figures  10  and  11  for  the  Lg/P  ratio  to  identify  earthquakes  as  such  most 
successfully  near  the  Gulf  of  Aqaba  and  least  near  the  Zagros.  Indeed  this  is  the  case,  as 
seen  in  Figure  12  where  we  have  plotted  contours  of  the  fraction  of  events  with  log(Lg/P) 
>  0  at  mb=3.5.  The  fractions  are  generally  high  where  the  observed  Lg/P  ratios  are  high, 
except  in  regions  where  the  mb(Lg)  threshold  is  high,  as  in  the  NE  part  of  the  Arabian 
peninsula.  The  fraction  of  identified  events  is  near  one  if  we  plot  them  at  mb=4.0,  as 
shown  in  Figure  13.  The  values  are  still  around  0.6  in  the  southern  Red  Sea,  due  to  poor 
propagation  and  local  values  of  log(Lg/P)  around  0. 


mejtcj.cdfl5.fAJ  .dIC.O.IdbyLgP 
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Figure  12.  Contours  of  fraction  of  events  identified  as  earthquakes  by  the  Lg/P  discriminant  at  mb=3.5  for 
the  proposed  IMS  network.  The  decision  line  is  log(Lg/P)  =  0.0. 
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Figure  13.  Contours  of  fraction  of  events  identified  as  earthquakes  by  the  Lg/P  discriminant  at  mb-4.0  for 
the  proposed  IMS  network.  The  decision  line  is  log(Lg/P)  =  0.0. 

We  do  not  currently  have  adequate  excitation  and  propagation  information  at 
higher  (>6  Hz)  to  accurately  simulate  the  regional  (Lg  or  Pg/Pn)  spectral  slope 
discriminants.  We  can,  however,  address  whether  higher  frequency  signals  could  be 
detectable  on  the  Arabian  shield  using  noise  spectra  obtained  from  series  preceding  signals 
in  Vernon’s  database.  These  stations  are  very  quiet  (Figure  5).  Using  the  propagation 
parameters  in  Table  2  and  a  standard  earthquake  source  model  (Brune,  1970),  we  find  that 
the  higher  frequency  signals  exceed  the  noise  out  to  1 0  Hz  (even  though  signal  power  is 
decaying  with  frequency,  the  noise  is  decaying  faster).  That  is,  At  locations  where  events 
can  be  detected  at  one  Hz,  they  can  generally  be  detected  at  ten  Hz.  If  the  station  planned 
for  Saudi  Arabia  is  comparable  to  Vernon’s  stations,  the  spectral  slope  discriminants  can 
be  applied  on  the  Arabian  shield. 

Each  of  the  stations  proposed  in  this  region  is  very  important  to  identification 
performance.  Due  to  lack  of  coverage  to  the  SE  of  SAUD  (in  the  Arabian  Sea),  this  station 
is  essential  for  monitoring  the  Arabian  shield  and  the  Sanai  peninsula.  Without  station 
LXEG,  there  are  no  regional  detections  in  NE  Africa  and  stations  THR  and  BEYT  are 
needed  for  regional  detections  south  of  the  Caspian  Sea. 
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2.  Maximum  Likelihood  Estimates  of  Teleseismic  mb  for  the  GSETT-3  Primary 
(Alpha)  Network 

2.1.  Introduction 

The  maximum  likelihood  estimates  of  magnitude  (Ringdal,  1976)  was  developed  to 
improve  measures  of  signal  amplitudes  when  are  below  detection  levels  in  some  or  all  of 
the  network  As  pointed  out  by  Von  Seggem  and  Rivers  (1978),  values  of  mb  derived  from 
arithmetic  mean  tend  to  have  a  positive  bias,  which  may  be  eliminated  by  the  use  of 
maximum  likelihood  estimates.  In  the  following,  we  simulate  a  suite  of  earthquakes 
recorded  at  the  GSETT-3  Alpha  network,  and  compare  the  mb’s  computed  using  the  two 
methods  with  true  values. 

2.2  Simulations 

The  simulations  described  here  use  the  Xnice  program,  described  in  detail  in 
Barker,  1996.  In  the  following,  we  specify  parameters  specific  to  this  problem  and  refer 
the  reader  to  that  document  for  methods,  models  and  other  parameters. 

2.2.1.  Source  Properties 

So  that  a  realistic  geographic  distribution  of  events  is  simulated,  we  extracted 
event  locations  from  the  Late  Events  Bulletin  for  1993  to  1996  from  the  prototype  IDC. 
This  included  about  20,000  events.  Seismic  moments  were  then  distribution  followed  the 
rule 


log  N„.  =  2.85-  log  M0, 
logM0>15, 

where  Ncum  is  the  cumulative  number  of  events  with  moment  >  M0  (note  that  in  Xnice  a 
moment  distribution  is  specified,  rather  than  a  magnitude).  The  stress  drop  were  governed 
by  log-normal  distribution  with  a  mean  of  10  MPa  (100  bars)  and  standard  deviation  of  5 
Mpa. 

2.2.2  Propagation  Parameters 

Amplitude  tables  are  from  the  Vieth  and  Claussen  tables.  Random  errors  in  log 
amplitude  of  teleseismic  P  due  to  propagation  followed  a  log-normal  distribution  with 
zero  mean  and  0.28  log-units  standard  deviation. 
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2.2.3  GSETT-3  Primary  Station  Network 


The  41  stations  used  in  this  simulation,  their  location  and  mean  log  noise 
amplitude  (microns)  at  1  Hz  are  given  the  following  table: 


Station 

Latitude 

Noise 

AKT 

50.434 

58.018 

-3.147 

ARAO 

69.535 

25.506 

-2.156 

ASAR 

-23.666 

133.905 

-2.774 

BDFB 

-15.633 

-48.000 

-2.899 

BGCA 

5.176 

18.424 

-3.349 

BJT 

40.018 

116.168 

-3.064 

BOSA 

-28.613 

25.416 

-2.726 

CMAR 

18.824 

98.947 

-2.920 

CPUP 

-26.331 

-57.329 

-2.919 

DBIC 

6.670 

-4.856 

-3.147 

ESDC 

39.675 

-3.965 

-2.751 

FIAO 

61.444 

26.079 

-2.488 

GECO 

48.836 

13.704 

-3.147 

HFSO 

60.134 

13.696 

-2.814 

KBZ 

43.950 

42.683 

-3.147 

LBNH 

44.240 

-71.930 

-2.375 

LOR 

47.268 

3.859 

-2.618  i 

LPAZ 

-16.288 

-68.131 

-3.552 

MAW 

-67.604 

62.871 

-2.303 

MBC 

76.242 

-119.360 

-3.347 

MIAR 

34.546 

-93.573 

-2.868 

MJAR 

36.542 

138.207 

-2.891 

NAO 

61.040 

11.215 

-2.450 

NPO 

64.771 

-146.886 

-2.007 

NRI 

69.400 

88.000 

-3.179 

PDAR 

42.780 

-109.560 

-3.315 

PDY 

59.600 

112.500 

-3.073 

PFCA 

33.610 

-116.460 

-3.005 

PLCA 

-40.731 

-70.550 

-2.955 

SCHO 

54.817 

-66.783 

-2.735 

SPITS 

78.180 

16.370 

-1.583 

STKA 

-31.882 

141.592 

-2.813 

TXAR 

29.334 

-103.667 

-3.243 

ULM 

50.250 

-95.875 

-3.015 

VNDA 

-77.519 

161.846 

-2.780 

WALA 

49.056 

-113.911 

-2.983 

WHY 

60.695 

-134.967 

-2.840 

WOOL 

-31.073 

121.678 

-2.919 

WRA 

-19.944 

134.341 

-3.020 

YKA 

62.493 

-114.605 

-3.433 

ZAL 

53.940 

84.805 

-3.463 
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The  GSETT-3  Primary  network  of  stations  has  a  wide  range  of  short-period  noise  levels 
ranging  from  a  high  of -1.6  (0.025  microns)  to  low  values  of -3.5  (0.003  microns). 

2.3  Results 

Figure  14  shows  the  locations  of  stations  and  events  used  in  the  simulations.  In 
Figure  15,  we  have  plotted  arithmetic  mean  mb’s  ( n%am )  versus  log  M0 .  Only  those 
events  which  meet  the  GSE  detection  criteria  are  included.  Also  plotted  in  this  figure  is 
the  ensemble  average  of  the  values  of  mb  within  0.125  of  the  plotted  log  M0  (the  averages 
are  indicated  by  x’s).  A  third  set  of  values  on  the  figure  are  the  values  of  mb  when  there  is 
no  propagation  variation  or  ground  noise  and  no  statistical  variation  in  the  stress  drop. 
We  refer  to  these  values  as  the  “true”  values.  The  curvature  is  due  to  the  change  in  comer 
frequency  with  moment.  We  see  that  mbam  exceed  the  true  values  by  about  0.5  at  the 
lowest  magnitudes  and  approach  the  true  value  at  higher  mb.  Figure  16  shows  a 
comparable  plot  for  maximum  likelihood  estimated  mb  ( mhmle ).  The  values  of  m"'k  are 
quite  close  to  the  true  values  across  the  magnitude  range.  n%am  and  mbmle  are  plotted 
against  each  other  in  Figure  17,  which  shows  that  the  rr^  am  exceed  the  mhmle ,  in  agreement 
with  Von  Seggem  and  Rivers  (1978). 


180‘ 120"  60'  0‘  60'  120‘180' 


180' 120'  60'  O'  60'  120' 1 80' 

Figure  14.  Locations  of  stations  (triangles)  and  events  (circles)  used  in  the  simulations. 
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Figure  1 5.  Arithmetic  mean  m,  (circles),  ensemble  average  m»(x)  and  true  m,  (diamonds)  are  plotted  versus 
log(moment). 
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Figure  16.  Maximum  likelihood  estimate  mb  (circles),  ensemble  average  mb(x)  and  true  mb  (diamonds)  are 
plotted  versus  log(moment). 
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Figure  17.  Arithmetic  mean  mb  are  plotted  against  maximum  likelihood  estimate  mb. 


The  incremental  number  of  events  versus  mhmle  and  mh  am  are  shown  in  Figure  18. 
In  each  case  the  values  are  normalized  by  the  total  number  of  events  detected.  The  peak 
number  of  events  for  mhmk  occurs  at  a  lower  value  than  1% am,  indicating  an  apparent 
lower  detection  threshold  for  mhmle .  The  cumulative  numbers  (number  with  values  greater 
than  mj,)  are  shown  in  Figure  19.  Also  shown  on  this  figure  are  the  true  values.  The  slope 
of  the  m"'le  distribution  at  high  magnitudes  (from  which  a  “b  value”  would  be  measured) 
has  about  the  same  value  as  the  slope  of  the  true  values,  while  the  n%  am  distribution  has  a 
greater  slope.  Note  that  the  slope  is  greater  than  one  (the  slope  of  the  log  moment 
distribution  above)  due  to  comer  frequency  effects. 
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Figure  19.  Cumulative  number  of  events  versus  mh"'le  (circles)  and  mh  "m  (diamonds)  and  true  mb  (inverted 
triangles).  Straight  dashed  lines  have  slopes  of -1  and  -2. 

These  results  indicate  that  replacement  of  mh  am  with  mhmle  as  the  network  mb 
values  in  future  bulletins  will  reduce  magnitude  bias  at  low  magnitudes  and  decrease  the 
slope  of  the  cumulative  distribution  of  events  versus  mb  in  accordance  with  conventional 
wisdom. 
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